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Abstract 

Based on a two dimensional odd-parity superconducting order parameter for 
Sr2Ru04 with p-w&ve symmetry, we investigate the single vortex and vortex 
lattice structure of the mixed phase near H&. Ginzburg-Landau calculations 
for a single vortex show a fourfold structure with an orientation depending on 
the microscopic Fermi surface properties. The corresponding extended Lon- 
don theory is developed to determine the vortex lattice structure and we find 
near H c \ a centered rectangular vortex lattice. As the field is increased from 
H c i this lattice continuously deforms until a square vortex lattice is achieved. 
In the centered rectangular phase the field distribution, as measurable through 
/iSR experiments, exhibits a characteristic two peak structure (similar to that 
predicted in high temperature and borocarbide superconductors). 
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I. INTRODUCTION 



With the observation of zero resistivity in Sr2Ru04 below T c = 0.93K, Maeno et a/.0 
discovered the first layered perovskite compound without Cu0 2 planes showing supercon- 
ductivity. Recent experimental and theoretical research provides considerable evidence that 
the superconducting state of Sr 2 Ru0 4 is of an odd-parity p-wave symmetry. 

Classifying the superconducting states of a generalized BCS theoryi, the mean fields 

A a> ,/(fe) = - Kw(fc,fc')( c fcv c -fcV') (1) 

k',cr,cr' 

are split into singlet and triplet pairing states making use of the antisymmetry condition 
A(fc) = — A T (— k). While singlet states are described through an even scalar function ip(k) 
in the form A(fe) = ia y ip(k), triplet states take the form 

A(fe) oc i(d(k) ■ &)& v , (2) 

where d(k) = —d(—k) is an odd vector functioni. Taking into account the crystal symmetry, 
the fc-dependent functions ip(k) and d(k) have to belong to an irreducible representation of 
the lattice point group, in our case the tetragonal point group D 4h . 

Initial proposals of odd parity superconductivity in Sr 2 Ru0 4 were founded on the itiner- 
ant ferromagnetism of related ruthenate compounds and on the similarities between Sr2Ru04 
and 3 He in the Fermi liquid corrections produced by electronic correlationsi. Recently, muon 
spin rotation (/iSR) measurement s@ have observed that a spontaneous magnetization begins 
to develop below T c . This finding is most naturally explained in terms of a time rever- 
sal (T) symmetry breaking state implying a multiple component representation of the 
point group. Furthermore, since the electronic structure of Sr 2 Ru04 is quasi-two dimen- 
sional, a gap function with a strong fc 2 -dependence is unlikely. Of the remaining two odd or 
even two-component representations of D^, the odd T^ u representation with basis functions 
{di(k), d 2 {k)} = {zk x , zk y } is the only without nodes at k z = and therefore is most likely 
to be realized in Sr 2 Ru04. 

In Ginzburg-Landau (GL) theory, we deal with order parameters f]j{R) depending only 
on the center-of-mass coordinate -R, defined through 

2 

A{R, k) = J2 r)j(R) i{dj{k) ■ a)a y . (3) 

i=i 

As usual we switch the transformation behavior from the basis functions dj(k) to the or- 
der parameter components rjj. The components (771,772) then share the rotation-inversion 
symmetry properties of (k x .k y ), and the broken T-state is (771,772) oc (1, ±z). Based on this 
order parameter, one of usS examined the vortex lattice structures of Sr 2 Ru04 near H c2 : 
for an applied finite field H along a high symmetry direction in the basal aft-plane, two 
vortex lattice phases have been found with a second superconducting transition between 
them, while a square vortex lattice has minimal energy when the field is along the c-axis. 
Band structure calculations@~i reveal that the density of states near the Fermi surface is 
mainly due to the four Ru 4d electrons, which hybridize strongly with the O 2p orbitals 
and give rise to three bands crossing the Fermi surface, labeled a, (3 and 7 (see Ref. 0). 
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Within the model of orbital dependent superconductivity!^, the orientation of the square 
lattice mentioned above depends on which of the Fermi surface sheets is responsible for the 
superconductivity. Recent experiments using SANS0 show a square vortex lattice aligned 
with the crystal lattice in a wide range of the phase diagram, indicating that the 7 sheet (due 
to the xy - Wannier functions) of the Fermi surface shows the superconducting transition!. 
Assuming that ferromagnetic spin fluctuations are responsible for the superconducting state, 
this finding is in agreement with 17 NMR experiments^! showing strong ferromagnetic spin 
enhancements only in the Ru 4d xy orbitals. 

In this paper we investigate the single vortex core structure and the geometry of the 
vortex lattice near the lower critical field H c \. We start from the Ginzburg-Landau theory 
for the two-component order parameter (771,772) and solve numerically for the vortex core 
structure of an isolated vortex. Next, we derive an extended London model valid for fields 
near if cl and large GL parameters k and determine the geometric structure of the vortex 
lattice. We find hexagonal lattices near H&, deforming into square lattices upon increasing 
the applied field. Taking the 7 sheet as relevant for superconductivity, the orientation of 
the square lattices is in agreement with experiments^. The individual vortex cores show 
square deformations. This vortex core structure is observable through scanning tunneling 
microscopy (STM) measurements since the in-plane coherence length of 900 A is well within 
the spatial resolution required by STM. Alternatively, Bitter decoration studies near if cl can 
be used to detect the predicted field dependent transition between hexagonal and square 
lattice, which would complete the picture found by SANS-measurements for higher applied 
fieldscJ. A further characteristic of the hex-to-square transition is the peak splitting in the 
field distribution calculated in Sec. [IV] and measurable in /xSR experiments. Note that a 
similar transition has been predicted in borocarbideEi'0 and high-T c superconductors^! and 
has been observed recently in the borocarbide materials^. 



II. GINZBURG-LANDAU RESULTS FOR THE VORTEX STRUCTURE 

We start from the Ginzburg-Landau free-energy density (/) for the two-dimensional 
representation T^ u of the tetragonal point group, with the Ginzburg-Landau coefficients 
determined within a weak-coupling approximation in the clean limit. These two approxi- 
mations seem reasonable for Sr2Ru04 given the ratios^ T c /eF ~ 10 -4 and l/£ ~ 8 (here 
T c denotes the critical temperature, ep the Fermi energy. I the mean free path and £ the 



coherence length). In usual dimensionless units, we hav< 



Ji8 



/ = ~(\v + \ 2 + \v-\ 2 ) + \ (\v + \ 4 + \v-\ 4 ) + 2 \v + \ 2 \v-\ 2 + \ i(v-v* + ) 2 + (v-v + ) 2 ] 

+B 2 + \Dr] + \ 2 + \D V _\ 2 
■ 1 + V [(D xV+ )(D xV _y + (D x r]S)(D x r] + )*} 

[(D yV+ )(D y ri„y + (D y r)„)(D y r) + y] 

[(D x rj + )(D y rj^y + (D yV+ )(D xV ^*} 

[(D xV+ y(D yV j + (D y rj + y(D x ^)}, (4) 
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where r)± = (rji ±ir} 2 ) /y/2, B = VAA, = —iV^/n — A^, f is in units B 2 /4ir, lengths are 
in units of the penetration depth A, B is in units \[2B C = $ /(27rA£), a = ao(T — T c ), and 
k = A/£. The anisotropy v = ({v*) — %(Vx v y)) I i( v t) + ( v x v y)) i s a parameter that measures 
the tetragonal distortion of the Fermi surface ((•) denotes averaging over the Fermi surface; 
note \u\ < 1 and v = for a cylindrical geometry). The free energy / is invariant under 
a simultaneous rotation of 45° around the c-axis and changing sign of v. This implies that 
all results derived below for v < can be transformed into the corresponding results for 
positive v > by a simple rotation, and we limit ourselves to the case v < 0. Further we deal 
only with external fields applied along the c-axis allowing us to omit all terms containing 
^-derivatives. 

At zero applied field the two degenerate solutions (r)+,r)-) = (1,0) and (r] + ,r]^) = (0, 1) 
minimize the free energy /. For nonzero applied field, the degeneracy is lifted, and depending 
on the direction of the applied field only one component r]± is stable. Throughout this paper 
we will fix the dominant order parameter to be r^_. The asymptotic vorticity of the dominant 
?7_-component then determines the field direction along the positive or negative z-axis. In 
the spatially inhomogeneous situations considered below, the second component r] + is driven 
through the mixed gradient terms [gradient terms such as (D x i] + )*(D y ri^)] and decays far 
away from the defect. Note, that to each solution B} a symmetric counterpart with 

reversed field direction exists 



{V-,V + ,B} - {<,<,-!?}. (5) 

Below, we first examine the asymptotic solutions of the GL-equations for v — 0. Second, 
we use the v = 0-solution and carry out perturbation theory in v to obtain the asymptotic 
solutions for v ^ (we ignore the coupling to the vector potential in the latter case). 

In the cylindrical case [y = 0) the problem exhibits rotational symmetry. Restricting to 
fields B = B(r)e z , the solution has the structure 

1-W = e»V(r), 
V+ (r) = e^ d V+ (r), 

A(r) = A(r)e e , (6) 

where (r, 6) are polar coordinates and n denotes the phase winding of the topological defect 
driving the vortex in the main superconducting component. The relative phase e l2d between 
r]_ and r] + corresponds to an angular momentum difference Al = 2 and comes in through 
the specific form (D x + iD y ) 2 of the mixed gradient terms for v — 0. The solutions for 
n — ±1 are the most stable ones. The small r expansion to cubic order takes the form 

n = l: ?7_(r) ~ [m\r - Qpg + ym? (l + ^-a^ 

A( r )~a?r-^i^r 3 , (7) 
8k 

n = -1 : rj-(r) - (m?r + y Qrf ~<+ ~/i {< + rf)) r 3 ) e "* 
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r 3 )e i 



V+ (r) ~ Mr + y (±m° - - ^a? (m° + p°)) 

AW^-^M!^, ' ' (8) 

where mj, p? 3 , and a? are the first expansion coefficients of r/_, 77+ , and A, to be determined 
by matching this small r expansion to the asymptotics at r — > 00. The latter reads 

^(r) = (1 -m(r))e me , 
7j + (r) = p(r) e l(n+2)e , 

A(r) = ±— + a(r), (9) 
«r 

where the ± stands for n — ±1. The asymptotic form of the GL-equations requires the 
three functions m(r), p(r), and a(r) to be of the form 



n + 2 



\2 N 



" +w ~ w e,< " +2,< '' (10) 

if we ignore the vector potential A(r) = 0. For the general case (A(r) 7^ 0) the main 
asymptotics is an exponential decay 

a(r), m(r),p(r) ~ e _r//A , (11) 

which is a consequence of transverse screening. From Fig. [l| we see, that for reasonably large 
k this exponent A approaches 1, which equals the usual penetration depth in dimensionless 
units. 

We proceed now with the anisotropic case v 7^ (here we neglect the vector potential). 
An expansion to first order in v produces the following result for small r 

n = 1 : 77_(r) ~ \ m\r ^ ^r 3 ! e ld + v f ^ — _ 2/ r e ' 

^W-^e^ + ^/e"' 9 , (12) 



n = —1 : V-{ r ) 



(rAr + g(p£ - 2m?)r 3 j e "* + < 3 r 3 e 3 ^, 



r? + (r) ~ (p°r + ^(m° - 2p°)r 3 ^) + z//° 3 r 3 e^ 9 , (13) 



and for large r 



V+(r) ~ e^+ 2 ) 9 + (e^+ 6 ) e + e^ e ) . (14) 
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The additional coefficients g± 3 and /° 3 are the first expansion coefficients of the higher 
angular momentum states (f± 3 from r/ + ~ e «(«+2±4)e anc j from ~ e «(n±4)0^ wri ich 
arise for z/ 7^ 0. Comparing the condensation energy for the two cases n = ±1 we find that 
the (n = — l)-case is the more stable: the admixed r] + order parameter rises linearly in the 
center rather than cubic as is the case for n — 1, and the asymptotic phase turns of r) + 
are the same in the center and for large r (which is also an argument against off-centered 
zeroes as seen for example in d^.^-superconductors). Both points lead to a higher gain in 
condensation energy, which should stabilize the (n = — resolution. Our numerical Ginzburg- 
Landau results (described in the next paragraph) for the critical fields H& = ein/Air confirm 
this finding: for k = 2.5 and v = —0.3, we obtain 

H cl (n = -1) = 0.435c, 
H cl (n = 1) = 0.47£? c , 

demonstrating the stability of the (n = — resolution compared to (n — 1). The correspond- 
ing s-wave result is H c i = 0A8B C . 

To obtain complete numerical solutions, we minimize the free energy density on a two- 
dimensional grid. The minimization is carried out by relaxing an initial configuration r] + = 0, 
?7_ = tanh(r) e ±td and A = iteratively. Indeed we find that the energetically favorable 
(n = — resolution shows no off-centered zeroes in A contour plot of the absolute value of 
the dominant (a) and ^-component (b) is displayed in Fig. |^. The admixed 77+ solution 
shows a fourfold structure with a maximal amplitude along the axes. In addition the vortex 
core center of the dominant order parameter is distorted with marked elongations along the 
diagonals of the crystal lattice for v < (and analogously along the axes for v > 0). This 
behavior is understood recalling the properties of the Fermi surface. For v < 0, the density 
of states (DOS) at the Fermi surface is larger along the axes than along the diagonals. A 
larger DOS however implies also a larger effective mass, which reduces the coherence length 
and thus the size of the vortex. Finally, the distribution of the absolute value of the _B 2 -field 
(see Fig. |3|) also shows clear square deformations aligned with the underlying crystal lattice. 

On the other hand, in the unstable (n = l)-case the admixed component r] + exhibits 
a negative phase turn around the center, and consequently four off-centered zeroes with 
positive phase turn are needed to match up this phase e~ l6 at small r to the asymptotic 
phase e 3ie at large r. This pattern is qualitatively well known from analog GL-calculations 
for s-admixture in d x 2_ y 2 high-T c superconductorsEl~i!. The dominant component r/_ and 
the 5 2 -field show similar features as in the stable case discussed above. 



III. EXTENDED LONDON THEORY 

We wish to determine the form of the vortex lattice close to H cl . Instead of solving 
the Ginzburg-Landau free energy directly, here we follow an idea recently developed in the 
context of d^.^-wave superconductors with a weak admixing of an s-componentM We first 
integrate out the admixed component, leading to a one-component free energy density. We 
assume k ^> 1 and !/< 1 and use London approximation to obtain an effective free energy 
density which depends only on the field B and which can be minimized with respect to the 
geometry of the vortex lattice. Though our results will be quantitatively inaccurate for the 
low-ft Sr 2 Ru04 material, we still expect to obtain a qualitatively correct picture. 
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The Ginzburg-Landau equation for i] + to first order in |^+|/|^-| reads 

1 - v 



= -77+ + 2 \ V . |V + D 2 n + + ^ (£>= - Dj) r/_ - — 



(D x D y + D y D x ) r]-. (15) 



Within London theory the dominant component f]_ has modulus unity (see Eq. (fL8|)) and the 
first three terms lead to the expression (1 + D 2 ) r) + . Solving formally for r/ + and substituting 
this expression into the free energy density, we obtain the effective one-component Ginzburg- 
Landau free energy 



f = ~\t]-\ 2 + -\v-\ 4 + \Dv- 



V(D,v)(l + D 2 ) 1 V(D, v) ?7_ + B 1 



with 



1 - v 
2i 



(D x D y + D y D x ) 



(16) 



(17) 



Here we have consistently neglected all terms of order r/ 4 , vrj 2 ^, and higher. 

Following the usual scheme for London approximation we replace the order parameter 
and its derivative by the superfluid velocity v 



V- 
Drj- 



[- V(j)(x) - A) j^ x) = -v e i4,{x) . 

K K 



(18) 



Keeping only terms up to second order in v and up to first order in the anisotropy parameter 
v, we obtain the London free energy in Fourier representation 



/(*) 



4k 4 l + k 2 /K 2 



(k 2 v 2 + 2u [(k 2 x - k 2 y ) (vl - v 2 ) - (2 k x k y ) (2 v x v y ) 



(19) 



Here, all quadratic expressions in v and B have to be understood as ViVj = Vi(k)vj(—k) 
(the results for vorticity n = ±1 are the same on this level of calculation; they would differ 
if one included higher orders in v). 

Next, we eliminate v. We obtain an equation connecting B and v by variation of the 
London free energy Eq. (|i~9D with respect to the vector potential A. The corresponding 
Euler-Lagrange equation reads 



iK 3 kAB = K 2 v-- 



4 1 + k 2 /K 2 



k 2 v + 2v (k 2 x - k 2 ) (v x x - v y y) - (2 k x k y ) (v x y + v y x) 



(20) 



We solve this equation for v up to first order in v and substitute the solution back into the 
London free energy density (Eq. (|T9"D), to obtain 

l + <? 2 f-, , 3 2 u ( (q 2 x -q 2 y) 2 (2q x q y ) 2 



f(q) 



1 

2 



l + q 2 K 2 - 



3 n 2\2 



1 + - q - - 
4 y 2 



B 1 



-- + J(q,v, k) B 2 , 



(21) 



where we introduced the short notation 

q = -k. (22) 

K 

The variation with respect to B(k) leads to the extended London equation J(q, u, K)B{k) = 
(in the limit K-^oowe recover the well-known result (1 + k 2 )B(k) = (no sources)). 



IV. VORTEX LATTICE CLOSE TO H cl 

We generate the vortex lattice by introducing the corresponding source terms at the 
vortex core positions into the London equations. A convenient form for these source terms 
isEH 

a(r) =Kj2e- K2{r ~ r > )2/2 , (23) 
j 

where the sum leads over all positions Vj of the vortex lattice. The magnetic field and the 
Gibbs free energy for the vortex lattice read 



B(r) = BJ2 



e iK.qjr e ~<lj/ 2 



where the specific form of the vortex lattice is given by the reciprocal lattice vectors qj. 
We parameterize the unit cell in real space through the angles 9\ and (9 2 , the ratio of the 
lengths of the unit vectors ( and the average field B (see Fig. |j), where the first three 
parameters describe the geometry and B the size of the cell. The lattice geometry depends 
on the applied field H and is determined by minimizing the Gibbs free energy Q{H) given 
in Eq. ( p4[) with respect to the unit cell parameters. Writing the operator J(q, u, k) in the 
form 



J(q, ,/,«) = ! + * 2 ( / + + /J 2)2 (V + f 9* - \ <7 4 cos40 g ) , (25) 



it is obvious that a change of sign in v still corresponds to a rotation in real and reciprocal 
space by 45°. We thus restrict the discussion again to negative values of v. The resulting 
vortex lattices exhibit a centered rectangular unit cell with the main axes aligned along 
the diagonals of the underlying crystal. Without loss of generality we thus assume Q\ = 
7r/4 — # 2 /2 and £ = 1 in the remainder of the paper. In Fig. [5] we plot the angle 82 
minimizing the Gibbs free energy for different values of v versus the field H (Fig. [|a: k = 
25.0, Fig. |5|b: k = 5.0). For large enough values of the anisotropy parameter u, the vortex 
lattice evolves from hexagonal at low applied fields to square at high fields. The lattices 
found in our calculations [y < 0) are always aligned with the underlying lattice (see Fig. ^|). 
For decreasing k, the crossover to a square lattice approaches to the lower critical field H cl 
and occurs at lower anisotropies z/, an effect which is believed to be even stronger if one 
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would include all nonlinearities into the theory. These results are in good agreement with 
recent experimental data of Riseman et al.0, who find an aligned square lattice in Sr 2 Ru04 
down to very small fields in the phase diagram. 

Apart form the vortex lattice symmetry, the shape of the individual vortices as well as 
the corresponding field distributions is of interest. The former, because STM-measurements 
are able to determine the form of the individual vortices directly, and the latter since the 
field distribution is the typical result of /iSR-measurements. While the shape of the I?-field 
around the vortex core is perfectly circular for v = over the whole range of the applied 
field H, we find square deformations of the 5-field for negative values of v (see Fig.^). This 
is in agreement with the single vortex results obtained in the previous section, which also 
exhibits the corresponding deformation of the fields around the vortex core. 

Experiments using ^uSR can observe the field distribution 



In Fig. [7|, characteristic results are shown for two = 25.0, H = 3.1 (a) and K = 

5.0, H = 0.5 (b). While at v — the field distribution shows only one peak with a small 
shoulder towards the minimal field value, the situation for v ^ is quite different. With 
increasing \v\ the v = 0-peak splits into two due to the appearance of two nonequivalent 
classes of saddle points in the B- field. Upon further increase of the low field peak 
gradually vanishes and finally develops into a broad shoulder. 



We have investigated the single vortex and vortex lattice structure of Sr 2 Ru04 near H c i, 
assuming a weak-coupling, clean limit description of the 2D odd-parity T^ u symmetry of the 
superconducting state in Sr 2 Ru04. The stable single vortex solution (denoted (1, —1)) takes 
the asymptotic form r]_) = (p(r)e t0 , [1 — m(r)]e~ ie ) . A second low energy vortex solution 
has the asymptotic form (?7 + ,^_) = (p(r)e 3ld , [1 — m(r)]e td ) (denoted (3,1)). Solutions of 
the latter type appear to be the minimal ones for o?-wave superconductors with an admixed 
s-wave component. Here, the symmetry allows for the competing (1, — l)-solution which has 
a lower energy as it contains fewer nodes and decays slower near the vortex core. The stable 
(1, — l)-solution shows clear four-fold deformations due to the anisotropy of the density of 
states at the Fermi surface. Current experimental results on the vortex lattice structure 
point towards an elongation of the vortex along the crystal lattice diagonals, in agreement 
with our theoretical analysis. A direct proof of such a structure should be observable through 
STM measurements. 

The structure of the vortex lattice has been determined using an extended London ap- 
proach where the admixed r] + component has been integrated out. The resulting London 
theory predicts that with increasing field the vortex lattice will continuously deform from 
a hexagonal to a square lattice. The orientation of the square vortex lattice depends on 
the anisotropy parameter v: For v < the square vortex lattice is aligned with the crystal 
lattice whereas for v > the orientation of the vortex lattice is rotated by 45°. Compar- 
ing to recent experimental results of Riseman et a/.0, who found the vortex and crystal 
lattices aligned in almost the entire field range, we conclude that v < 0. The continuous 




(26) 



V. CONCLUSIONS 
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hex-to-square transition should be observable through Bitter decoration experiments. An 
alternative signature of this transition is given through the splitting in the peak of the field 
distribution P(B). These field distributions should be observable in /iSR experiments. 
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FIGURES 



FIG. 1. Exponent A governing the decay of v = solutions. 

FIG. 2. Contour plot of GL-calculations for the absolute values of the dominant t/_- (a) and 
the admixed r/ + -component (b) for the parameters n = 2.5 and v = —0.3. The contours are 0.99, 
0.975, ... for (a) and 0.03, 0.045, . . ., 0.225 for (b) 

FIG. 3. Contour plot of the absolute value of the |-B|-field for the same parameters as in Fig. ^. 
Notice the strong square deformations of the vortex core along the diagonals of the crystal lattice. 
The contours are 0.03, 0.06, . . ., 0.42. 

FIG. 4. Schematic illustration of the real unit cell. 

FIG. 5. Vortex lattice structure as function of the applied field H. The GL-parameter is k = 2.5 
(a) and k = 25.0 (b), and 02 denotes the angle between the unit cell vectors minimizing the free 
energy density given in Eq. (p4[). 

FIG. 6. Contour plot of the field distribution for v = — 0.4, k = 2.5, H = 0.25. The horizontal 
and vertical axes of the plot correspond to the axes of the underlying crystal lattice. The field 
distribution shows an aligned square lattice with weak square shaped deformations of the vortex 
cores. 

FIG. 7. Field distribution P(B) for different values of the anisotropy v at k = 2.5, H = 0.3 (a) 
and k = 25.0, H = 3.1 (b). 
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